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AN ITERATIVE GUIDANCE SCHEME AND ITS 


APPLICATION TO LUNAR LANDING 


SUMMARY 


QG&8 

A guidance scheme for vehicle flight from lunar orbit to a prescribed point on a 
spherical, non- rotating moon is presented. The equations of motion have been simpli- 
fied only to permit a closed solution for the thrust magnitude and thrust direction. The 
trajectory computations themselves are made under more realistic and accurate assump\ 
tions and are not included. 

This scheme is another approach to the problem of adaptive guidance mechaniza- 
tion for vacuum flight. The more salient features of the scheme are the limited number 
of presettings, the closeness to optimization, the closed form solution, and the ’’homing” 
feature; i.e. , the guidance information is improved as the end point is approached. 


The effectiveness of the scheme and the required thrust variations are displayed 
in Tables I and n. Two flight profiles have been considered: continuous burn of one 
RL-iO engine for main braking from a 100 km lunar orbit to the beginning of the termi- 
nal braking phase at 300 meters altitude; ignition of three RL-10 engines at the 15 km 
periselenum of a Hohmann transfer ellipse from a 100 n. mi. circular orbit with burn to 
300 meters altitude. A description of a timer ignition scheme is given in the appendix. 


SECTION I. INTRODUCTION 



It has been shown by Fried [1] and Lawden [2] that the optimum thrust direction 
in a vacuum over a flat surface with a constant gravitational field, constant thrust, and 
constant specific impulse is 


9R/ay T + (T - t) 3R/9y T 

tanx = 8R/8x t + (T - t) 3R/3x t 


a + bt 
c + dt 


(1) 


where a, b, c, and d are constants depending on the boundary conditions. For example, 
if the velocity vector alone is constrained, the direction of this vector can be achieved by 


tan x = 


8R/8y T 

9R/3x t 


= a’ 


(2) 





V m 1 

tt = -V ; • = - T, 
m ex’ m 

where irij represents instantaneous mass. 

If only the end velocity vector is specified, it has been shown that a constant 
thrust attitude direction, x > will achieve t-his end condition in a near-optimum fashion. 
Equations (2. i) can then be integrated over T seconds in closed form: 


x(T) = x t + (cos x) V ex ln^Y 

(2.3) 

y (T) = y t + (sin x) V ex In — ~ - gT. 

Solution of equations (2.3) for x yields 


Smt) -Ji.*rr . (2 . 4) 

x(T) — x, 

Inspection of Figure 1 shows the same result geometrically. 

The magnitude of the velocity vector is obtained through the burning time T; i. e. , 



AV. - *J(£(T) -x 1 ) 2 +[J( T) -yj + gT] 2 , (2.5) 

which may be set equal to the characteristic velocity 

AV. = V In — • . (2.6) 

l ex t - T ' ’ 


Equations (2. 4) , (2. 5) , and (2. 6) can be solved for and x for any combination 
of current state velocities x lt y t and required end velocities x(T) , J^(T). Since current 
state variables are changing, the T-computation and x-determination proceed stepwise 
using new Xj, y t values as they are obtained. 

So far, only an end velocity vector condition has been obtained. In order to 
constrain an end altitude also, it is necessary to introduce a linear thrust attitude law, 

X = x - K, + Kjt , 
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which has been justified through experience, is used in developing the scheme. These 
laws make it possible to generate the function in flight using as information state vari- 
ables of velocity and displacement, thrust, mass flow, and desired end conditions. 

The scheme is modified to approximate a spherical moon assumption by intro- 
ducing an average gravity direction and magnitude between the instantaneous state and 
the desired end condition. No longer are a’ , a t and bj constant throughout the flight as 
in the unperturbed flat moon case but constant over the length of time between evalua- 
tions. Indeed, one might think of a series of flat moons between the instantaneous point 
and the end point", an assumption which becomes more valid as the end point is ap- 
proached. These modifications to approximate a spherical moon do not destroy the 
ability to solve for a* , a ls bj, and T in closed form. 

% 

Range constraint is achieved by the throttling capability of the engines. The 
required thrust is inversely proportional to the required range. The change in thrust 
affects, ;of course, the solution of x and T. These quantities are updated, however, to 
be compatible with the new thrust. 


SECTION H. DESCRIPTION 


The principles of the scheme can be displayed most clearly by assuming a two- 
dimensional flat surface as a model. In a subsequent section, the extension to a homo- 
geneous spherical body, still in two dimensions, is presented. Finally the specific 
application to a lunar landing mode is shown. 

' The equations of motion relative to a flat surface in a constant gravitational field 

are 


x = F/m cos x 

( 2. i) 

y = F/m sin x - g, 

where F/m represents acceleration due to a constant thrust and a constant specific 
impulse, and x represents the direction of the acceleration vector against the horizontal. 
Mathematically expressed, 

^ = £•— . ( 2 . 2 ) 

m m mj + mt ex r - t 

Since vacuum flight and constant specific impulse are assumed, F and m are 
constant; m < 0. The following definitions also prevail: 
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where 


G = y t + y,T - sin x V gx |j - T In 7 - 7 ^ - ^ - y(T) (2. 12) 

S = V ex cos H T, -f ln< rM • 

SECTION III. EXTENSION OF THE GUIDANCE SCHEME TO A 
SPHERICAL BODY ASSUMPTION 


The theory can be extended to the spherical body assumption by assuming an 
average value for the gravity direction 0 * and magnitude g* , where 

- 0i 

0* = — 7 — ; (3.1) 

gi + g T 

g* = — g • (3.2) 


Figure 2 depicts the coordinate system used in the guidance computation. The 
£-17 system, with origin at the desired injection point, is formed by rotating the x-y 
system , space-fixed at the launch site through the terminal range angle <p^: 


4 = x cos 0 T - (R q + y) sin 0 T 


( 3. 3a) 

17 = x sin <P T + ( R q + y) cos . 

A similar relation holds for the velocity components. 

4 - x cos 0 T - y sin 0 T 


V = x sin (j> T + y cos <P T . 


( 3. 3b) 
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FIGURE 2. FLIGHT GEOMETRY 
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The steps taken in deriving the thrust attitude angle are the same as those 
displayed in Section II. 


A i = l T " il " g* T sin <£* 
A n = n T - n i + g* T cos </>*. 


* -1 M. 

y = tan 1 — . 

h ai 


AV. = V In 


i ex t - T 


(3. 4) 


(3. 5) 


(3.6) 



(3. 7) 


X| = X^-K, + K 2 t (3.8) 

x = X - K 1 + K 2 t- <*> T . (3.9) 


The equations of motion in the 1-17 coordinate system then becomes 
l = F/m cos x^ + g* sin 0* 
n = F/m sin x^ - g* cos 0*. 

Again the coefficients Kj and K 2 appearing in (3. 8) are determined to preserve 
the velocity condition imposed by equation ( 3 . 5) while enforcing a condition on altitude. 
Thus, 

?(T) = r, 1 - g*T cos <t>* + sin ? V ex In (3. 10) 
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^(T) = T? t - g*T cos <£* + sin x^ V gx In - ^ T 

- K > 006 * I V ex ln 7^F + K ’ V ex cos * { [‘ T + T ln TTY] 
77(T) - f(T) = 0 


t?(T) 


= 17 T = *?i + V lT - g* y cos <p* + V ex sin x ^ |t - (r - T) In 
-K,V ex cosT { Jt-(t-T) 

+ V ex K J co S T 4 |-|-r[(T-T)ln 7 i T - T ]} . 


Solving equations ( 3.12) and ( 3. 13) simultaneously for K t and K 2 , 


Kj = — J 


G 

S 


K 2 = 


G . r 
— In 


S r - T 


where 


J = r In — - T 

r - T 


A 


G - vt + iiT - sin x { V ex Tj - T In - g* y °° e ** - l(T) 


> 


S = V cos x 
ex A £ 


[ TJ -T to <TTY»] • 


(3.11) 

(3. 12) 

(3.13) 

(3. 14) 

(3.15) 

(3. 16) 
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where X , Kj, and K 2 preserve the velocity condition already obtained, while also en- 
forcing a desired altitude condition. Equations (2. i) become 


x = F/m cos ( x - K t + K 2 t) 
y = F/m sin (x - Kj + K 2 t) - g. 


(2.7) 


hi order to preserve the velocity condition, Kj and K 2 must be chosen so that the 
y obtained from the equations (2. 7) and y are the same. It is assumed that 
cosf-Kj + K 2 t) « 1 and sin(-K 1 + K 2 t) » -Kj + K 2 t. 


y = yj + ( sin x - K t c os x ) V gx In - ■ + K* cos x V fix [-T + t In 


ex 


F r + Tln 7^rJ _gT - 


( 2 . 8 ) 


y( T) - y = 0 = -K t cos x V ex In — - + K 2 cos x 


ex 




(2.9) 


The coefficients Kj and K 2 must be chosen also to satisfy the altitude condition 


y(T) : 


“ — 

y(T) = y t + y t T + (sin x - K 2 cos x) V gx T - (t - T) In 


+ K 2 cos x V 


ex 


ln 7T 


- T 


gT 2 

2 * 


( 2 . 10 ) 


Equations (2. 9) and (2. 10) can be solved simultaneously for K t and K 2 : 


K. =f J 


K ’ = f ln <7TT» 
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It may be helpful to go into some detail on the computation of T, the time-to-go, 
which is determined from equations ( 3 . 6 ) and ( 3. 7) , since their form would seem to 
indicate the need for an iterative procedure to solve them. However, an iteration 
computation in flight is undesirable since there is no guarantee of convergence in the 
allotted time. Therefore, the following method is employed: Let T' be a guess of the 
time-to-go so that 


T = T' + <5T. 

Then 

A £*' = i T - li - g* T ’ sin 0* 

At?' = ?7 t “ Vi + g^T' cos <P* » 

and 

av ! 2 = \_h ‘ " g * T ' sin + 

- fi i + g*T- cos 4 > 2 . 

Now, referring to equations ( 3. 4) , (3. 7) and ( 3. 17) , 

AV 2 = [ 4 " *1 ‘ g * (T ’ +6T) Skl ^*] 2 + 

[^ T - + g*(T’ + 6 T) cos 0*] 2 , 

so that 

AV? - AV ^ 2 = 2A^’ g* cos <f>* - 2A| ' g* sin + g* 2 (<5T) 2 . 
For convenience of expression, the following relation is introduced: 
A = ( A 17 T g* cos (p* - A£ ' g* sin </>*) , 


(3.17) 


(3. 18) 


(3.19) 


(3. 20) 


( 3 . 21 ) 


( 3 . 22 ) 
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bo that 


AV* * g* 2 ( 6T) 2 + 2X6 T + AVV* . 


(3. 23) 


The total final inertial velocity in the guessed time T’ is 


=V 4 . +V Q In 


Tj li ex T _ t ’ 


or 


(3. 24) 


AV. = V In 


. — V 1U ___ • 

i ex t - T 


(3. 25) 


Expressing (3.25) as a series approximation. 


AV. - - V ex In (l - —) - V + i 

i ex \ t J ex T - T' 2 \T - Ty 


(3.26) 


Let 


L = V In 


ex t - T* 

V, 


K =— ; 

■t - T’ 


then equation ( 3. 26) becomes 


AV t * L + K6T + ~K • — » L + K6T, 

A ^ T - 1 1 


(3. 27) 


if second order terms are discarded. Squaring equation. ( 3. 27) and setting it equal to 
equation ( 3. 23) gives 


g* 2 ( 6T) 2 + 2X6T + AV' 2 = L 2 + 2KL6T + K 2 ( 6T) 2 . 


( 3. 28) 


ii 


Collecting coefficients of like terms results in 


- a( ST) 2 + 2bST + c = 0 


(3.29) 


where 


a = (K 2 - g* 2 ) 

b = A - KL (3.30) 

c = AV'. 2 - L 2 . 

l 

Solving the quadratic for ST , 


ST 


b + »Vb 2 + ac 
a 


(3.31) 


The positive root is chosen in equation (3.31) , since b is negative, and ST should be 
small. The time to go, T, may now be determined and the correct velocity components 
become 


a£ = A|* - g* ST sin <p* 
Aft = A^' + g*6T cos 0*. 


( 3. 32) 


The spherical Earth assumption assumes a knowledge of the terminal range angle 
( prj , . Presetting a range angle based on previous experience is one approach; however, an 
on-board computation of (p^ is more desirable since it will change accordingly with 
perturbations which may be encountered. Therefore, the following approach is taken. 

The horizontal distance obtained from a linear acceleration on a flat surface in 
time T* is 

X T = f ( y i + v ex ln TTt") dt ’ (3, 33) 

o 
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or 


x T = xi + T* 


h 


, T . „ T - T* , T - T*" 

+ v„„ + V ov - ^r~ In — 


ex ex t* 


(3. 34) 


Dividing by the radius vector r,p. 


^rp = + 


r T L 


„ „ „ T - T* , T - T* 

V< + V + V In 

1 ex ex T* t 


(3.35) 


Referring to equation ( 3. 1) , 




(3.36) 


This approximation has proved quite sufficient and has the characteristic of 
improving in accuracy as the end condition is approached. 

Briefly, in retrospect, the values preset are final altitude, h-p; final velocity 
components, £-j. and r^p and physical constants of the central body; values obtained in 
flight are instantaneous state variables x lf y lt x t , y 1# and some measurement of the 
specific impulse or fuel flow rate. 


SECTION IV. APPLICATION OF THE GUIDANCE SCHEME 
TO LUNAR LANDING 


The guidance equations developed in the previous sections for ascent are virtually 
identical to those needed for lunar descent. The two essential modifications are a change 
in coordinate system due to the way the instantaneous state information is obtained and 
the use of range prediction to vary thrust, since it is desired to arrive at a prescribed 
spot on or above the lunar surface. 

The potential presenoe of a beacon lends itself to the coordinate systems shown 
in Figure 3. The £-17 system has its origin at the beginning of the terminal braking 
phase; the £-axis passes through the instantaneous vehicle position (M) . The beacon 
provides slant range, D, and slant range rate, £>. It is also assumed that the angle 0 
and its angular rate 0 are known or can be computed. 
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* The beacon is located at the origin of a space-fixed system with its y-axis 

parallel to the radius vector to that point. The coordinates of the predesignated terminal 
braking point and the vehicle velocity at that point in the space-fixed system are 
(Xiji» y T> Xijij y*ji) . 

The vector extending from the terminal braking point to the vehicle is 


| = i (D cos 0 - x^,) + j (D sin 9 - y^.) 


(4.1) 


| = i cos 0p + j sin 0p. 


(4. ia) 


Using these relationships, 

— — i , , i 

£i = (4 • 4) 2 = K D cos 0 - x^.) 1 + (D sin 0 - y T ) z ]2 , 

D cos 0 - x T 
COS 0p = 7 , 

li 


sin 0 p = 


D sin 0 - y T 


tan 0p = 


D sin 0 - y 


D cos 0 - 


*T 


(4. 2) 


(4. 3) 


(4.4) 


The velocity components in the D-0 system are transferred into the space-fixed 
system : 


Xj = D cos 0 + D 0 sin 0 

(4.5) 

y t - - f) sin 0 + D0 cos 0. 

The instantaneous values are obtained by transferring to the £-17 coordinate system: 

L = x. cos 0p - y< sin 0p 

(4-6) 

t/ 1 = Xj sin 0p + yj cos 0p. 
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Similarly, for the terminal velocity components, 


£ t = *T COS 0p " Sin 0P 
V T = x T sin 0p + y T cos 0p. 


(4. 7) 


A significant change necessary is the determination of the average gravity direc- 
tion 4 >* , since the 77-axis is not parallel to the local vertical at the terminal braking 
point. 


<P 


* 



+ 4> t ) - op. 


(4. 8) 


Equation ( 4. 8) can be more readily understood by considering first the 77-axis parallel 
to the local vertical at the terminal braking point which yields equation ( 3. 1) and then 
subtracting the angle ( 0p + $ T ) necessary to rotate the local vertical into the 17- axis. 

The guidance equations are identical to those developed in the previous two 
sections. The form of equations (3.14) and (3. 15) is the same; however, fhe terms 771 
and tj(T) which appeared in equations (3. 16) are dropped, since in the £*-77 coordinate 
system for lunar landing, they are zero. Thus , 


J = 


r In 


T 

T - T 


- T 


(4.9) 


G = ^iT - \ g* cos <t > * T 2 - V ex sin x [j - T In 77^] 

S=V cosx T J ~-7T In (--- ■ ) . 

ex A l_ 2 t-T J 

Equation (3.9) becomes the following to relate the thrust altitude angle to the 
space-fixed x-y coordinate system : 


X ~ X~ Kj + K 2 t + 0p. 


(4.10) 
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•Thrust Control 


In the previous sections , no constraint has been placed on range. For lunar 
landing, however, it is desired to control range since a specific point on the lunar 
surface should be reached. To accomplish this task, variable thrust is superimposed 
on the guidance scheme. 

The thrust ( F) is inversely proportional to the distance , the predicted distance 
corresponding to the thrust ( F) . 


F = , (4.11) 

where K is a proportionality constant. There will exist a desired thrust level ( F*) such 
that the required distance ( £) is covered during the remaining time. 



(4.12) 


Solving yields 



Equation (4. 13) is the thrust control law for selecting the thrust level. The displacement 
£* is the predicted distance that will be traversed during the remaining flight at the 
thrust level (F) : 


£* = !l-|(T). (4.14) 

The total average acceleration in the ^-direction is 

*|(t) = — cos x + E* sin ft*. (4.15) 

The value of cos x is an effective average value of cos x^ over the time T. The 
positive sign on the gravitation component is due to the negative ft*. 


17 


Using the fundamental laws of calculus, 

T T 

/ d$(t) = |(T) - = / 4<t) dt 


$(T) = |j+ / |(t) dt 


t t 

/ d|(t) = |(t) - It = / |(t) dt 


|(T) = |t + / *|(t) dt. 


Combining equations (4. 15) , (4. 16) , and (4. 17) and placing into (4. 14) , 

T f r t f f 1 1 

$* = -/ ) ii + f mi + ^t cos * + g * sin ** dt [ dt - 


Performing the designated integration yields 


k = - 4l T - V p X cos x T - (t - T) In 


— ]- 
T - T J 


* . .* 

g* sm<p — 


Using the defining equation ( 3. 5) for \ - X and the defining equation ( 3. 6) for 
AV, the relation 


AV = - At— = y i n - T 

cos X ex t - T 


is obtained. 


Placing equation (4.20) into equation (4.19) yields 


I* -= - |iT - V ex cos x T + (t - T) A| - g* sin (p* -j 


(4.21) 


Using 


A| = krr "il - g* T sin 4>* 


gives 


• • T 

= tA£ - T ( - g* T sin </>*) - V ex cos x T - g* sin (p* — 


Rewriting equation ( 2. 6) and substituting in equation ( 4. 20) yields 

a£_ 


T = r \ 1 - e 


V ex cos X 


(4. 22) 


(4. 23) 


(4. 24) 


For convenience, let 

.. .A£ 


z = 


V cos x 
' ex A 


(4. 25) 


For computational purposes it is worthwhile not to have the exponential in ( 4. 24) ; 
therefore, equation (4.24) is expressed in a Maclaurin series: 


Z 2 Z 3 Z 4 Z s 


T- t - t I 1 - Z + - - ^r + 4T- 5 , + ••• 


or 


■•fM HI}- 

Placing ( 4. 27) into ( 4. 23) , and using ( 4. 25) , gives 




(4. 26) 


(4.27) 


• 1 * 

- - — g 

»T 2 & 


1 

2 


( 

r- 

** 

. z 

. z f . z\ 


i _ — 

1 3 



l 


0 


(4.28) 
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Hence, the value of £* needed to select the proper thrust magnitude using equations 
(4.13) is obtained. 

The task remains to determine cos x so that equation (4. 25) may be evaluated. 
Let cos x be defined as follows: 


2 T t 

cos X = f f cos dt 2 , 
o o s 


(4. 29) 


where x^ = x^ - K i + K 2 t. Performing the designated integration, 


2 i i 

cos x = p / p [sin x - sin (x - K t ) ] dt, 
o 


(4. 30) 


cos x = l ~ ^ [ cos x - cos ( x - Kj) ] - T sin (x - Kj) l . 


(4.31) 


Expanding the cos x term and rearranging, 


c °s x = ^Fp (cos (x-K t ) [1 - cos K 2 T] - sin (x - Kj) [K 2 T - sin K 2 T ] 


(4. 32) 


From the computer viewpoint it is most desirable to eliminate as many trigo- 
nometric functions as possible. Although cos x can immediately be obtained from 
equation ( 4. 32) the following steps are taken to obtain a more suitable form for the 
computer. Using the trigonometric identity, 

cos x = A cos ( x - Kj + oo) = A cos ( x - Kj) cos co - A sin ( x* - Kj) sin co. ( 4. 33) 


Equating terms with 4. 32, 


A c°s co = [1 - cos K 2 T]. 


A sin co = [K 2 T - sin K 2 T] . 


(4. 34) 
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Squaring ( 4. 34) , adding, and expanding into series , we obtain 


A 2 = 4./ r (K |T >! _ ( k 2T)‘ 

K 2 4 T 4 j [2 24 

a 2 _ _ 4 _ 1 

K 2 4 T 4 1 L 4 48 

A =n/i-(|-|) (K 2 T) 2 




(4. 35) 


Using the binomial expansion gives 


A = 


_ (K 2 T) 2 
36 


(4. 36) 


and expanding ( 4. 34) in a series gives 


& = sin' 


2 

A K 2 2 T 2 


(K 2 T ) 3 

6 



(4. 37) 


A further simplification can be made if it is recalled that the enforcement of 
equation ( 3. 12) implies approximately that 


/ 


o 


T 

Xdt * / (x-K t + K 2 t) dt 
o 


XT * (X- K, + K 2 ^)'T. 


(4. 38) 


Thus, 

K 2 T « 2Kj . (4.39) 

Figure 4 shows this approximate relationship. 
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X - Kj + K 2 T 



x - Ki 

L 

t seconds 

FIGURE 4. THE APPROXIMATE RELATIONSHIP OF K 2 T « 2K t 


Evaluating ( 4. 33) , using ( 4. 36) , ( 4. 37) , and ( 4. 39) , yields 


K< 


Ki 


cos x = l 1 - — J cos U - 7 


3 ’ 


(4. 40) 


With these relationships , the required thrust can be computed to fulfill the range con- 
straint. It should be noted that changing the thrust magnitude affects the solution of 
K l5 K 2 , x » and T; however, these quantities are updated with the new thrust evaluation. 
The following block diagram may be helpful in showing the input and the mechanization 
of the scheme itself. 
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x, K„ K 2 , T 


x “ X T } F* which affects control quantities above 



FIGURE 5. BLOCK DIAGRAM SHOWING MECHANIZATION OF SCHEME 
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SECTION V. DISCUSSION OF RESULTS 


Two rather typical flight profiles for lunar landing to exemplify the effectiveness 
of the guidance scheme are ( 1) continuous burn of one RL-10 engine, producing a 
nominal thrust of 15,000 lbs, from a 100 kilometer circular lunar orbit to a terminal 
braking point of 300 meters altitude, and ( 2) main braking by three RL-10 engines, 
producing a nominal thrust of 42,000 lbs, from the 15 kilometer periselenum of a 
Hohmann transfer ellipse to a terminal braking point of 300 meters altitude. 

Tables I and II show the absolute value of the thrust variation and the displace- 
ment and velocity dispersions evaluated at the 300 meter terminal braking point for 
various initial dispersions in state variables and engine performance, and for radar 
errors emanating from the beacon. The guidance parameters were determined every 
ten seconds. The trajectories generated by the scheme adhere closely to those gener- 
ated by variational calculus procedures. The maximum payload difference between the 
two methods is approximately 50 lbs. 

Since it is believed that the choice of initial errors is conservative, it can be 
seen that the scheme itself has virtually no error and is as accurate as the source of 
information it requires. Indeed, the scheme can be made as accurate as desired by 
increasing the number of evaluations limited, however, as the desired end point is 
reached, by the possibility of control parameters becoming indeterminate. 

A crude analysis has been made of computer requirements to implement the 
scheme. These requirements are listed in detail in Table III. It should be borne in mind 
that the computer time estimates are based only on the program used in checking out the 
scheme. It may be of interest to note that the guidance evaluation takes ten seconds for 
one step on the RPC 4000, which means, for the cases in Tables I and II, that the co- 
efficients would be obtained in real time. On the GE-225 digital computer the evaluation 
can be made in 0. 5 seconds. The scheme requires 1300 words of storage; however, it 
is estimated that with optimum programming, this requirement can be reduced by 
approximately 50 percent. 

The single RL-10 case is somewhat unrealistic in that 15,000 lb nominal thrust 
has been assumed, and thrust increases for some perturbations cannot be realized with 
the existing engine. The case does show, when compared with the 3-engine case, that 
the percent thrust variation required of the single engine is considerably higher than 
that required for the 3-engine case for nearly identical disturbances. In view of this 
trend, and because of long powered flight times and ranges required of the single engine 
flight profile, the single engine descent cannot be highly recommended. 
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TABLE I 


SCHEME AND RADAR ERROR ANALYSIS 


Continuous Burn of One Engine from 100 km Orbit 


W t = 71,000 lb 

f nom = 15 ’ 000 1b 

W T * 43,655 lb 


Perturbation 

|AF| 

lb 

Ax 

m 

Ah 

m 


Ay 

m/s 

AWgt 

-lb 

0 

736 

.05 

-. 04 

-. 001 

-. 001 

0 

Ax = 25 km 

709 

.03 

-. 02 

-. 002 

-.002 

-287 

Ay = 25 km 

1064 

.07 

-.03 

.000 

-.010 

-606 

Ax = -50 m/s 

730 

.07 

-. 05 

.001 

-. 013 

-194 

Ay = 50 m/ s 

791 

. 06 

-.03 

.001 

-.011 

-277 

Algp — —4. 25 

1826 

2. 57 

2. 54 

.254 

.566 

-144 

AWgt = 500 lb 

1704 

3.21 

2. 74 

.285 

.640 

-155 

Radar Errors 











1% Slant Range 

1288 

-2. 86 

-5. 37 

-.001 

-.015 

184 

1% Range Rate + 
1% Slant Range 

502 

-.29 

-2. 22 

.101 

.259 

6 

10% Angular Rate 
of Slant Range 

1185 

10.81 

-5.51 

1. 036 

-1. 253 

265 

1° A Slant Range 
Angle 

1372 

-1. 08 

-3. 93 

.414 

- .595 

233 

RSS of Scheme 
Error + Radar 
Error 

3850 

12 

10 

1.2 

1.7 

875 
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TABLE II 


DESCENT FROM 15 km PERISELENUM OF 
HOHMANN TRANSFER ELLIPSE WITH 3 ENGINES 


F NOM 


= 42,000 lb 


Wj = 71,000 lb 
W T = 46,543 lb 


Perturbation 

1 af | 
lb 

Ax 

m 

Ah 

m 

Ax 
m/ s 


AWgt 

-lb 

0 

64 

-.01 

.01 

-.013 


0 

Ax = -25 km 

193 

O 

O 

• 

.00 

-.020 

-.009 

52 

Ay = 10 km 

196 

02 

.03 

-.028 

-.013 

-391 

Ax = - 50 m/ s 

59 

O 

o 

• 

-.01 

-.018 

-.012 

-523 

Ay = 50 m/ s 


-.02 

.03 



-194 

AI sp = -4.25 

Wm 

13. 66 

5.43 

■sa 

.666 

-31 

AWgt = 500 lb 

2194 

13. 84 

5. 48 

1. 530 

.690 

-31 

Radar Errors 











1% Slant Range 

3139 

-7. 76 

-4.91 

-.016 

.002 

12 

1% Range Rate + 
1% Slant Range 

1895 

10.51 

4. 38 

1.317 

.603 

-28 

10% Angular Rate 
of Slant Range 

689 

18. 78 

-43. 76 

1. 105 

-2. 726 

-6 

1°A Slant Range 
Angle 

1457 

o 

CO 

• 

1 

-17.60 

.311 

-2. 320 

9 

RSS Scheme 
+ Radar 

5200 

30 

48 

2.8 

3.8 

690 
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TABLE IB 


COMPUTER REQUIREMENTS FOR GUIDANCE SCHEME 

L Built in addition, subtraction, multiplication and division. 

IL Subroutines 

A. , sin 

B. cos 

C. nT 

D. log 

E. tan -1 

F. integration for gravity computer, if required 
ITT, Storage Requirements 

A. integration 16 (erasable) 

B. storage and program 910 

C. subroutines 250 

TOTAL * 1300 

IV. Time Requirements per Step, no printing 

A. GE 225, internal floating point 0. 5 sec 

B. RPC 4000, programmed floating point 10. 0 sec 
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APPENDIX 


IGNITION CRITERION FOR MAIN BRAKING 


It is possible that some cases of main braking may require ignition before the 
beacon has been acquired. For this reason, .the timer approach is attractive both for 
simplicity and accuracy. The following assumptions are made: 

1. The period of the actual lunar orbit in which the vehicle flies prior to 
braking is known - measured from the point at which the slant range rate, D, is zero. 

2. The ignition time for the nominal case, t ig , is known and is available aboard 
the vehicle. 


3. The thrust attitude program and the thrust magnitude for the nominal main 
braking case are precalculated and are available aboard the vehicle. 

The time of ignition for any case, tj , is then defined as 


ig 


- P - P + 


t ig 


(a) 


At this time the vehicle begins to fly, with the nominal thrust attitude program and the 
nominal thrust magnitude, until the beacon is acquired, when the guidance scheme pro- 
duces the necessary thrust vector to accomplish the mission. 

The following table displays a typical example, based on a nominal main braking 
phase from the 20 kilometer periselenum of a Hohmann transfer ellipse, the aposelenum 
of which is at 100 kilometers. It is assumed that three RL-10 engines with a nominal 
thrust of 40, 000 pounds are available and that acquisition of the beacon occurs sixty 
seconds after ignition. 

The table shows the change in thrust magnitude and direction called for when the 
beacon is acquired due to errors in the original ellipse, and errors resulting from fly- 
ing the nominal thrust magnitude and direction in the region where the beacon is not 
available. The errors in state variables at this point can easily be handled by the 
guidance scheme. 
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ERROR ANALYSIS OF PERTURBED ORBITS 
FOR TIMED IGNITION 


Nominal Thrust = 40, 000 lb 

Guidance Scheme Becomes Effective 60 Seconds After Ignition 


km 

hp 

km 


AF60 

lb 

Ax 

km 

Ay 

km 



•AX60 

AWgt 

lb 

100 

20 


0 

0 

0 

0 

0 

0 

0 


10 


135 

-.98 

-10. 23 

6. 79 

3.11 

13. 02 

-59 


30 


642 

1.41 

10. 26 

-6. 73 

-3. 47 

-13.51 

-66 


20 


-91 

.28 

.04 

-2. 20 

-.55 

.10 

0 

110 

20 


73 

-. 29 

-. 04 

2.18 

.56 

-.11 

0 
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